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way, the scalar propagator of the usual effective theory. 
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1 Introduction 

The usual way to approach scalar quantum electrodynamics (SQED) is by per- 
forming the electromagnetic minimal coupling in the free Lagrangian of Klein- 
Gordon (KG) scalar field theory Q. An alternative way is to start from the free 
Duffin-Kemmer-Petiau (DKP) Lagrangian instead of the KG one. 

The free DKP theory is a theory for scalar and vector fields [|2|, |^, [|, |5], |B| given 
by a Lagrangian formally similar to that of spinorial QED. The fundamental 
differences are the algebraic relations satisfied by the (3^ matrices in DKP theory, 
which play the role analogous to 7 M in spinor QED. 

It is known that in the free field case the DKP and KG theories are equivalent, 
both in classical and quantum pictures H> H> 0- However, there are still no 
general proofs of equivalence between these theories when interactions and decays 
of unstable particles are taken into account 0, [TTJ] (in this context see also 



references [|12| 13], which relate different results by using both DKP and KG 
formalisms with strong interactions). Some progress in this direction has been 
made recently. For instance, it was shown that both theories are equivalent in 
the classical level for the cases of minimal interaction with electromagnetic || [14| 
and gravitational J0j] fields . Strict proofs of equivalence between both theories 
were also given for the cases of interaction of the quantized scalar field with 
classical and quantized electromagnetic, Yang-Mills and external gravitational 
fields [0,0. 

Perhaps one of the most evident advantages in working with this theory is the 
fact that derivative couplings do not appear between DKP and the gauge field 
(this property has been used by Gribov recently, who employed the vector sector 
of DKP theory to study the quark confinement problem |17|]). Such property will 



result in manifestly covariant expressions for the interaction Hamiltonian and the 



vacuum expectation values of time ordered products of fields. Another advantages 
are the formal similarity with spinor QED (which facilitates the adaptation to 
the scalar case of technics developed formerly in the spinor case [|T2| , |T8|| ) and 
the fact that this formalism allows an unified treatment of the scalar and vector 
fields. 

One of the difficulties in working with SQED based on KG equation (SQED- 
KG) is the presence of a term of second order in the coupling constant in the 



interaction Hamiltonian, which causes trouble in proving gauge invariance fl9 |. 
In SQED based on DKP theory it was achieved, by an effective approach, that this 
second order term does not contribute to S-matrix, and thus it can be neglected 



when we construct the Feynman rules for the theory || [16| . 

In this paper we will consider SQED-DKP in the framework of Epstein-Glaser 
causal perturbative method. This method was formulated to give a mathemati- 
cal rigorous treatment of ultraviolet divergences in quantum field theory. In this 
framework such divergences do not appear anywhere in the calculations due to 
the correct splitting of causal distributions into its advanced and retarded parts 
[ f20| , |2i]| . Our goals are to obtain a non-effective and mathematically well defined 
theory for SQED-DKP and recover the results of the corresponding effective the- 
ory already obtained in the usual perturbative approach. In addition, this work 
must be viewed as an initial step in the attempts to rigorously establish the 
renormalizability of the theory which is still an open question. 

This paper is organized as follows. In section 2 we present the DKP theory. 
In section 3 we review the main features of the Epstein-Glaser causal method 
and present the basis to construct the second order S matrix. In section 4 we 
consider Compton scattering and address the question of how the propagator of 
the effective theory emerges in the causal approach. In this context the gauge 
invariance plays a crucial role. In section 5 we calculate the scalar vacuum polar- 
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ization tensor. In section 6 we calculate the self-energy and, by using an Ward 
identity, also the vertex correction in the limit of zero momentum transfer. In 
section 7 we make our concluding remarks. The implications of gauge invariance 
are presented in the Appendix. 



2 Duffin-Kemmer-Petiau theory 

The free DKP theory is given by the Lagrangian P, ^ [| 

C=lW^rf-mW>, (1) 
where ip is a multicomponent wave function, ip = ip^r/ , and rf = 2 ((3 ) 2 — 1. (3 P 

are a set of matrices (/i = 0, 1, 2, 3) satisfying the algebraic relations 

Z^/?"/?' + = ^g vp + (3 p g^. (2) 

The equations of motion are then 

{%^d v -m)ip = and ^ (i(3 p d v + m) = 0. (3) 



It is known [2~2~ that the algebra (Q) has only three irreducible representations, 
whose degrees are 1, 5 and 10. The first one is trivial, having no physical content. 
The second and the third ones correspond, respectively, to the scalar and vectorial 
representations. In this work we shall restrict us to the scalar case. 

Defining ^= and using relations (Q) it can be shown that, for any four- 
vector q the following relation is satisfied 

f(f - q 2 ) = 0. (4) 

By using this relation and the plane wave solutions for the equations of motion 
we can verify that p 2 = m 2 . 



As an ilustrative example of an explicit representation of @ we can get 
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Using this specific representation we can construct explicitly projection operators 
and show that not all components of the field ip are independent. For a detailed 
explanation we refer to references |J |8| . Here we only quote the explicit form of 
ip in this representation: 



/ 



(□ + m 2 )^ = 0. 



(6) 



\ V^V ) 

Thus we readily see that not all components of ip can be independent. Actually, 
we can choose only two of these components as independent ones. We shall make 
use of this explicit representation in the remainder of this work only for ilustrating 
the arguments, being our results valid for any scalar representation of relations 



Now we apply the standard procedure of canonical quantization in the free 
Lagrangian ([!]) and obtain || 

bPa ( x )^ + b(y)} = \ s ab( x - v) and \}Fa{ x )^b(v)\ = -] s ba(y - x ) > 

where ip~ and ip + contains only annihilation and creation operators, respectively 
and 

S±(x) = ^m + m)UA±(x), (7) 



A ± (x) = / d 4 p 5{p 2 - m 2 )6{±p°)e- ip - x (8) 



where 

A±< - X) = (ar). 

are the positive (negative) frequency parts of the Pauli- Jordan distribution 

A(x) = A + (x) + A~(x). 
As it is well known, this later distribution has causal support ETJ] and can be 



written as 

A(ar) = A rct (x) - A adv (x), (9) 

where A ret (x) and A adv (x) have, respectively, retarded and advanced supports 
with respect to the point x. Analogously we define 

S(x) = S + {x) + S~{x). (10) 

We can see directly that this distribution also has a causal support, since it is 
a linear combination of derivatives of A(x), and the differentiation of a causal 
distribution does not affect the causal property of its support. Now, by ([?[) and 
@, it is possible to write 

S(x) = S rct {x) -S adv (x) , (11) 

where 

S rct (x) = — [i3(i3 + m)]A vct (x) and S adv (x) = — [i3(i3 + m)} A adv (x) . (12) 
m m 
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As we will see in the next section, the above splitting of S(x) in retarded and 
advanced parts is not the unique possible. 

The interaction with electromagnetic field is introduced by the minimal sub- 
stitution — > — ieA^ in the Lagrangian ([!]), which becomes 

C = C F + Ci; (13) 

£f = % ^X^-^- (14) 

A = e^3"V4,, ( 15 ) 

where we have used e > 0. 



3 The Epstein-Glaser causal approach 

In the Epstein-Glaser's causal method [20 the ^-matrix is constructed with- 
out any reference to Hamiltonian formalism, its explicit form being obtained by 
making use of certain physical conditions - with causality playing a major role. 
In this approach the S-matrix is viewed as an operator-valued distribution given 
by the perturbative series 



°° 1 f 

S(g) = 1 + 2J — / d A xi . . . <fx n T n 

n=l U ' J 



(xi,...,x n )g(x 1 )...g(x n ), (16) 



where g(x) is a c- number test function supposed to belong to the Schwartz space, 
g(x) G 5(R 4 ). The symmetric n-point functions T n (X) (X = {xi, . . . ,x n }) are 
the basic building blocks to be inductively constructed, from the knowledge of 
Ti(x), by means of the requirements of causality 



S(gi + g2) = S(g 1 )S(g 2 ) , 



if supp gi > supp g 2 , 



(17) 
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and translational invariance 

U(a, l)S(g)U{a, l)' 1 = S (g(x - a)) . (18) 

In the above equations U(a, 1) is an usual representation of the Poincare group 
V\_ in the Fock space and the notation supp g\ > supp g2 signify that all points 
in the support of #2(2) occur at times previous than all points in the support of 
9i(x). 

In terms of the n-point functions T n these causality and translation invariance 
conditions are given by 



if {^l; • • • j •^m} ^* \^m+l) ■ ■ ■ > ^n} > 



and 



U(a, l)T n (xi, . . . ,x n )U(a, 1) 1 = T n (x 1 + a, . . . , x n + a) , (20) 

respectively. 

Making use of these requirements we are able to construct the T n , order by 
order, from the explicit form of T\(x). It is known that, based on general argu- 
ments such as correspondence JJJ , Ti = iCj , where C^p is the term of first order 
in coupling constant in the interaction Lagrangian and is written in terms of free 
fields H . 

Now, let us sketch the inductive procedure (for a detailed account see |^T], p3|j). 



Suppose that all T m (X), with m < n — 1 are known, then one can define the 
distributions 



P 2 



(21) 

Pi 

where P<i stands for all partitions Pi : {x±, . . . , = X U Y, X ^ into 



disjoint subsets with | X \= m, \ Y \< n — 2. In (21) T n (X) refers to the n-point 
distributions corresponding to a series for the S'~ 1 -matrix analogous to ( Jl6| ) and 
can be obtained by formal inversion of S(g), giving 

n 

UX) = J2(-YY, T ^ X J ■ ■ ■ T -™' ( 22 ) 

r=l P r 

where P r indicates all partitions of X into r disjoint subsets: X — X\ U . . . U X r , 
Xj ^ 0, \Xj\ = rij. Of course, since all T m (X), m < n — 1 are given by the 
induction hypothesis, also the T m (X) with m < n — 1 are known. 

If in (^j]) the sums are extended in order to include the empty set X = we 

get 



pO 
1 2 



(23) 



-^nv^lj ' " " ' "^w) "I - • • • > In)? 

where P 2 ° stands for all partitions P° : {xi, . . . ,0^-1} = X U Y. A glance at 
equations (|2"3"| ) shows that A n and i? n are not known because they contain the 
unknown T n . However, the distribution defined by 
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D n {x u ...,x n ) = R n -A n = R' n - A' n , (24) 

is known. 

Making use of causality it turns out that R n has retarded support and A n has 
advanced support, i.e. 

supp,R n (X) C r+_ 1 (a; fl ), suppA n (X) C T'^Xn), (25) 

with 

Fn-i( x ) = {(xi,...,x n _i) | xj e T^(x),Vj = l,...,n- 1}, 

(26) 

^(x) = {y\(y- x f > 0, ±(i/ - x°) > 0}. 

The distribution D n has causal support, supp-D„ C T^-^UT'^. In fact, a general 
proof of the causal support of D n only exists for n > 3 - for n = 2 we must verify 
this explicitly. Then, decomposing D n in advanced and retarded distributions we 
obtain the T n distribution by fl23"|) . 

The operator- valued distributions which we shall have to split are of the form 

Xij . . . , x n 

k j I m 

where ip, if) are the free boson fields of DKP theory and A stands for the free 
gauge boson fields. In this expression d\ are numerical tempered distributions, 
d\ £ S (R 4n ), with causal support. Because of the translation invariance, it is 
sufficient to put x n = and consider 
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d{x) = d k n (x 1 ,...,x n . 1 ,0) E S'(R m ), m = 4n-4. (28) 

The nontrivial step is the splitting of the numerical causal distribution d in 
the advanced and retarded distributions a and r, respectively. From the fact 
that r + (0) PI r~(0) = {0} we can see that the behaviour of d(x) in x= (or, in 
momentum space, p = oo) is crucial in the splitting problem. With this in mind, 
a classification of the distributions is given in which d{x) e S (R m ) is called 
singular of order u if its Fourier transform d{p) has a quasi-asymptotics do(p) at 
p = oo 



21, 23 



lmp(6)(dg),4>(p)) = (d ,4>), (29) 

V 

(4> is the inverse Fourier transform of 0), with power-counting function p(5) sat- 
isfying 

for each c > 0. Of course, there is an equivalent definition in the coordinate space 



|PH|, but, since the splitting is more easily performed in the momentum space, 



this one is sufficient for our purposes. 



Then, we have two distinct cases 20, 21 : i) uj < - in this case the solution 



of the splitting problem is unique and the retarded distribution can be defined by 
multiplication by step functions; ii) u > - now the solution can be no longer 
obtained by multiplying d by step functions and, after a careful mathematical 
treatment, it may be shown that the retarded distribution is given by the "central 



splitting solution" pl | 



i f r d(tp) , s 

Hp) = — / dt- r \P- r. 31 
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This solution has the very important feature that it preserves the symmetries 
of the theory, in special Lorentz covariance and gauge invariance. However, in 
contrast with the case uo < 0, the solution of the splitting problem ([H]) is not the 
unique one and, in momentum space, the general solution is given by 

rip) = rip) + Yl C ^ ( 32 ) 

[a[=0 

where the C a are constant coefficients which are not fixed by the causal structure 
- we need additional physical conditions in order to determine them. 

In expression (^) use is made of the minimal distribution splitting condition 
which says that the singular order cannot be raised in the splitting. This condi- 
tion, crucial for a correct prediction of the anomalous magnetic moment in QED4 
21|1 and in the analysis of the dynamic mass generation in (2+1) dimensions 



24], |2l)|, will also be useful here. 



Now we apply the inductive steps above to construct the two-point distribu- 
tions for SQED-DKP theory. Then, the one-point distribution for the DKP field 
interacting with electromagnetic field is given formally by % times the interaction 
term flTSP in the Lagrangian (|T3|), 

T x {x) = xe : ^(x)P^(x) : A^x) = -f^x) , (33) 

where all fields entering in this expression are free fields and e is the physical 
charge. The normal ordering is necessary in order to have a well defined expression 
for the product of field operators at the same point. 

To go from n = 1 to n = 2 we take (^) and construct the distributions 
A' 2 (xi,x 2 ) and /2' 2 (xi,a; 2 ) from expressions (|23|). So, we have 

A' 2 (x 1 ,x 2 ) = e 2 : ^(x 1 )^(x 1 ) :: ^(x 2 )^(x 2 ) : A^ Xl )A v (x 2 ); (34) 
ti 2 {x u x 2 ) = e 2 : ^(x 2 )^(x 2 ) :: ^(x 1 )f3^(x 1 ) : A^ Xl )A u (x 2 ) , (35) 
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and then 

D 2 (x l ,x 2 ) = e 2 {:^(x 2 )(3^(x 2 ) ::^(x 1 )(3^(x 1 ) : (36) 
- : ^{x 1 )(3^{x 1 ) :: ^(x 2 )^(x 2 ) ^Ap{x x )A v {x 2 ). (37) 

By using Wick's Theorem this expression can be written as a sum of terms, 
each of them consisting of a product of field operators and Wick contractions 
2lj| . In this work we shall be concerned to the terms corresponding to Compton 



scattering, vacuum polarization and self-energy. From these we will determine 
the vertex correction via an Ward identity. These terms are given by 

D^(x u x 2 ) = e 2 ^ cd 

x{: ip^xjipdfa) :: A^x^A^x^ : ^ B (x 2 )ih{xi) - ipb{xi)ip c (x 2 )) 

+ : il) b {xi)il) c {x 2 ) :: A lx {x 1 )A u {x 2 ) : [^(^2)^(^1) -^(^1)^2)]} ; (38) 

D^( Xl ,x 2 ) = e 2 (3^ cd :A,( Xl )A u (x 2 ): 

x[^ c (x 2 )i) h {xi) i) d {x 2 )i) a {xi) -ip a (x 1 )ip d (x 2 )] , (39) 

Df\x u x 2 ) = e 2 ^ cd 

x{: 'ip a (x 1 )tp d (x 2 ) : [^ e (x2)if) b (xi)A v {x 2 )A (t (xi) 

- ^6(^1)^0(^2) A^(xi)A u (x 2 )} 

+ : ^6(^1)^0(^2) : [^d(x 2 )ip a (xi) A v (x 2 )A^(xi) 

-^ a {x x )^ d {x 2 ) aJ^£(x 2 )}} , (40) 

where the Wick contractions are defined as 

^M$ b {v) = [VC(z),^(v)] = 7<s£(*-v); ( 41 ) 
Mx)Mv) = \^{x)M{y)] = -~^ a {y-x) ; (42) 
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A^AM = [M x )' A t{v\ =>ig»vD+(x-y) , (43) 

and S^ b (x) and S ba (x) are given by (0). Dq{x) is the positive frequency part of 
the zero mass Pauli- Jordan distribution, Dq{x) = -^p- / d 4 p 5(p 2 )9(p°)e~ ip ' x . 



4 The Compton scattering 

We denote respectively by D\(xi,x 2 ) and D^(xi,X2) the first and second terms 
inside curl brackets in fl3"E|). By using fl4"T| ) and E^) we have 

D\(x l ,x 2 ) = ie 2 (3^ b (3 v cd : tfj a (x 1 )tfj d {x 2 ) :: A M (xi) A,(x 2 ) : 

x{5S;(x 1 -x 2 )+5+(x 1 -x 2 )} (44) 
= ie 2 f3^ b f3 v cd : V , a( aJ i)V , d(^2) :: ^(^1)^(^2) : S& c (xi - x 2 ) (45) 

where we used the definition fllPl). The first term inside curl brackets in (f44|) 
cames from B%(xi,X2), whereas the second cames from ^(xi, a^)- Similarly, 

L> n (xi,x 2 ) = ie 2 l3^ h (3 u cd : ^ 6 (xi)V> c (£ 2 ) " ^(x^A,^) : 

x {-5+(x 2 - x x ) - S da (x 2 - Xl )} (46) 

i^PabPcd ■ ^b{xi)ip c {x 2 ) :: A M (xi)A„(x 2 ) : 
*{S da {x 2 -xi)}, (47) 



where again the first term inside curl brackets in ([46]) cames from R' 2 ll (xi,x 2 ) 
and so on. As we said in Section 2, the distribution S(x) has causal support. 
Then, we have verified explicitly that the distributions ( [I5| ) and (57) have causal 
support too. 

Since S(x) itself is the numerical distribution we must split into retarded and 
advanced parts in equations (ff5]) and fl4"7|), a splitting solution is trivially obtained 
from dTT|) . But, since S"(x) has singular order u; = -as we can verify by (|29|) and 
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( |30D - this splitting is not unique. So, accordingly with (j32|) , the general solution 
for the retarded distribution in configuration space is 

r( x ) = S mt {x) + C5(x), 

where C is an arbitrary constant. Now we construct the numerical distribution 
^(ari,^) from @, 



t 1 (x 1 , x 2 ) = r(xi,x 2 ) - r'(x 1} x 2 ) , (48) 
where r'(x\,x 2 ) is the numerical distribution associated with R^(xi,X2), which is 



obtained from the first term inside curl brackets in (f44|) . Then, 



t 1 (xi,x 2 ) = S rct (xi - x 2 ) - S (xi - x 2 ) + C5(xi - x 2 ) 

= -5 F (x 1 -x 2 ) + C5(xi-x 2 ) , (49) 

where we have defined 

- S F (x) = S TCt (x) - S~(x) = S adv (x) + S + (x) = -—iana + m) A F (x) , (50) 

m 

where A F (x) is the usual Feynman scalar propagator. 

In a similar way we find (note that the reference point for splitting is x 2 ) 

t l \x u x 2 ) = S adv (x 2 -x 1 ) + S + (x 2 -x 1 ) + C'5(x 2 -x 1 ) 

= -S F (x 2 -x l ) + C'5(x 2 -x 1 ) . (51) 

The constants C and C will be determined later by the requirements of charge 
conjugation and gauge invariance. 

Now the complete two-point distribution for Compton scattering is given by 

T 2 Compton (xi, x 2 ) = T 2 \x!,x 2 ) + T 2 n (xi, x 2 ), (52) 
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where 

T 2 I (x 1) x 2 ) = ie 2 : ^(x 1 )(3H\x u x 2 )(3 v il;(x 2 ) :: A ll (x 1 )A v (x 2 ) : 

and 

T 2 II (x 1 ,x 2 ) = *e 2 : ^{x 2 )(3 u t 1I {x 1 , x 2 )(3^{x l ) :: ^(xO^fo) : . 

It is very simple to require charge conjugation invariance of the theory [|| |2l| . 
Here we only quote the condition that arise from this invariance, namely 

C = C. 

Thus we can immediately see that 

Tl{x u x 2 ) = T 2 n (x 2 ,xi), 

what implies that T 2 ompton (xi, x 2 ) is symmetric in its arguments, as it would be. 

In the Appendix we calculate the remaining constant C by requiring gauge 
invariance of the theory. There we find (see eq. (|9C|)) 

m 



where I is the 5x5 identity matrix. Turning this result into ( |49| ) and ([51]) we 
obtain 

t l (xi,x 2 ) = -{S F (x 1 — x 2 ) 5(xi — x 2 )\ 

m 

= t l \x 2 ,X!) . 

We denote the distribution inside the curl brackets in the above expressions as 

T c (x) = S F (x) - -5(x) . (53) 
m 

It is direct to see that this distribution is the Green function for the DKP equation, 
i. e., 

{i$-m)T c (x) = 5(x) . (54) 
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So, in this sense the distribution (p3|) is the correct "propagator" for the DKP 
scalar particle, that is the same as the "effective" propagator of references || [16 
For later use, we write this propagator in the momentum space 



T c (p) 



1 



(27r) 2 m 



j>{$> + m) 



(55) 



p2 _ m 2 _|_ 

With the above results we write finally the two-point distribution for the 
Compton scattering: 



T? omptoa ( Xl ,x 2 ) 



-ie 2 {: ^(xi)^T c (x! - x 2 )(3^(x 2 ) : 



+ : i>{x 2 )(3 u T c {x 2 - x 1 )^(x 1 ) :} : A^x x )A v {x 2 ) : (56) 



5 The vacuum polarization 

Now we consider the term (|39D, associated with scalar vacuum polarization. After 
substituting the explicit form of Wick contractions, this term is given by 

D^ ac {x 1 ,x 2 ) = e 2 : A„{ Xl )A v {x 2 ) : 

xTr{^S~( Xl - x 2 )(3 v S + {x 2 - Xl ) - [3 u S~{x 2 - x l )(3^S + {x l - x 2 )} , 

(57) 

where the first term inside curl brackets cames from R'(x\ — x 2 ) and the second 
from A'(x\ — x 2 ). Making y — X\ — x 2 and defining 

P^{y) 4 -e 2 Tr{^S + (y)(3 v S-(-y)} , 

we write fl57|) in the form 

D^ c (x u x 2 ) = {P^(y) - P*"H/)} : M*i)M**) : • ( 58 ) 

This distribution has a causal support, as required. To see this we use the fact 
that S(x) has a causal support and write the term inside curl brackets in the 
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above expression as 

^[S{vWS + {-y) - S + {-yWS{-y)\ . 

Before splitting the causal distribution above, we consider the Fourier trans- 
form of P' MU (y), 



= ^Tr J dy^S + (y)^S-(-y)e 



ik.y 



Substituting into this expression the explicit forms of S + (y) and S (—y) given 
by (|7|) and (|5p, we have 

(59) 

x6(p°)5(p 2 - m 2 )6(k° - p°)5[(k - p) 2 - m 2 ]. 
To compute the trace inside this integral we need the following trace properties 

Tr^ 1 /^ 2 .../^ 2 - 1 } = ; 
with n = 1, 2, 3, .... So, we obtain 



Tr{/3^ + m)/3 v (^-i>)(^-^-m)} = m 2 [Ap IM p v - 2(k^p u + p^k") + k^k u ] . 

When substituting this result into (H^), and taking into account the distribution 
S[(k — p) 2 — m 2 ], we observe that P^ u {k) is a second rank symmetric tensor which 
satisfies 

k,P^{k) = , 
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what means that the vacuum polarization term (58) is gauge invariant. Using 
this fact, we can write, as usual, 

P^(k) = {k»k v - k 2 g^)B(k 2 ) , 

where 

B{k2) = w^ {k) ■ 

The computation of P^{k) is straightforward and the final result for P^ u (k) is 

= wff" r ) ~ 1 6 {k2 ~ 4 '" 2) 6ik0) ' (60) 

Now, in momentum space the numerical distribution associated with fl58"D , 
which we must split, is given by 

—e 2 I ' k^k v \ ~ 

P^(k) - P^(-k) = ^ ^_ - g^j d{k) , (61) 

where 

nk 2 ( 4m 2 



d(k) = — \1 - -p-J (k 2 - Am 2 ) sgn(fc°) . (62) 

As the tensor character of (|6i"D does not matter in its splitting procedure, the 
problem reduces to the splitting of d(k). To do this, we first determine the 
singular order u> of this distribution by using ( p9|) and (p0[). So, we find u = 2. 



To simplify the calculations we make use of the fact that k in ( |62| ) is timelike. 
So, we can choose a Lorentz frame such that k = (k°, 0). In this case the central 
splitting formula Q3l|) can be written as plf : 



(p°-i0) 3 (k°-p° + i0) 
Substituting (|62"D into this integral and making the substitution s = (p ) 2 , we 
obtain 

i ,. u I /" . 1 ( _ 4m 2 \ 2 



f (fc°) = - (k°Y { P.V. / ds^^- 2 1 - 

4m 2 
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-27TS, 



gn(fc>[(A; ) 2 -4m 



4m 



9 \ 2 
2 \ 2 



(A; ) 2 V (A; ) 2 

The distribution r'(k°) comes from the term —P^i—k) in fl58|). So, 



f'(A; ) 

6 V (A; ) 

Now the two-point distribution is given, in an arbitrary Lorentz frame, as 



vr(A; ) 2 ( 4m. H 



t(k) = f{k)-f\k) 



12 J s{k 2 -s + iO) 

Am 2 



1 - 



4m 



2X1 



Finally, the two-point distribution for the vacuum polarization in configuration 
space is given by 



T 2 Vac (x 1 ,x 2 ) = -i : A^{x 1 )W u (x 1 -x 2 )A v (x 2 ) : 



where 



W v (k) 



—ie 2 ( k^k 



(2tt) 4 V k 2 
1 /W 



(2vr) 4 V A; 2 



<T I t(fc) 
«r) ft(Ar) 



and 



n(jfc) = —A; 4 
v ; 12 



ds 



Am 2 



s(k 2 - s + iO) 



4m 2 \ 2 



(63) 



For < k 2 < Am 2 , the result of this integral can be writen as 



fl(k) 



—e 



4m 2 -A; 2 n . 1 
-3p-(ft»M-l) + 5 



J& and 9 6 (0, f ; 



where sin 2 ^ 

For spacelike k we get 



IL(k) 



e 2 m 2 



12A(1-A) 



;i + A) 3 lnA + ^(l-A 3 ) 
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where k 2 = — ^—^-m 2 , and < A < 1. This result can be analytically extended 
for k 2 > 4m 2 by making — 1 < A < and taking InA =ln|A| — in. 

These results are not the most general solutions for the splitting problem, 
as we had to split a distribution with singular order uo = 2. The most general 
solution tl(k) is given by 

U(k) = tl(k) + C + + C 2 k 2 , 

where the normalization constants Co, C M and C 2 are not determined by causality. 
They are determined from the requirements of a zero mass for the gauge field, 
parity invariance and the identification of m with the physical observable mass. 
The procedure is the standard one and we simply quote the results: Co = C M = 
C 2 = 0. 

6 The self-energy and vertex corrections 

The self-energy of DKP particle is given by (|40l) . After substituting the Wick 
contractions we obtain the following two terms: 

Df i {xi,x 2 ) = -e 2 :^(x 1 )^[S-(x 1 -x 2 )D+(x 2 -x 1 ) 

+S + ( Xl - x 2 ) D+ ( Xl - x 2 ) (3^(x 2 ) :] ; (64) 

D!f{x u x 2 ) = e 2 : ^{x 2 )^ [S + {x 2 - x x ) D+ {x 2 - Xl ) 

+S~ (x 2 — Xi) Dq (xi — x 2 ) Pfj,ip(xi) :] . (65) 

Therefore we note that 

D^f(x 1 ,x 2 ) = -Df { (x 2 ,x l ) . (66) 

Let us first consider Df elf (xi,x 2 ). From it we can determine the two-point 
distribution Tj self (xi, x 2 ) . From 

A sclf (^i^ 2 ) = E! 1 (x 1 ,x 2 )-A[(x 1 ,x 2 ) 
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= Ri(x 1 ,x 2 ) - Ai(x!,x 2 ) 

and 

Djf (xi, x 2 ) = R'n fri, x 2 ) - A' u (xi, x 2 ) 
= Ru(x 1 ,x 2 ) - A u (x 1 ,x 2 ) , 

and using (|66|) we conclude that 

R' n (xx,x 2 ) = A[ (x 2 ,xi) ; 
A' n {x 1: x 2 ) = R[(x 2 ,xi) ; 
Rn(xi,x 2 ) = A I (x 2 ,x 1 ) ; 
An(xi,x 2 ) = Rj(x 2 ,xi) . 

So, we have 

T I sclf (x 1 ,x 2 )=T I f f (x 2 ,a; 1 ) (67) 

and conclude that the two-point distribution for the self energy T self (xi,x 2 ) = 
Tj self (x\, x 2 ) +Tfj (xi,x 2 ) is symmetric, as it would be. Moreover it suffices to 
compute one of the distributions, say Tj self (xi, x 2 ), and the other will be given by 

(0)- 

Returning to Df 1 ^ (xi, x 2 ) in (|64D, we must first verify the causal property of 
the support of this distribution. The term into brackets can be writen as 

S (xi - x 2 ) (x 2 - xi) + S + (xi - x 2 ) D (xi - x 2 ) 

and, as S (xi — x 2 ) and D (xi — x 2 ) are causal distributions, so it is also 
Df { { Xl ,x 2 ). 

The numerical distribution we have to split is 

d x {y) = -e 2 ^ [S-(y)D+(-y) + S + (y)D+(y)] ^ . (68) 
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We denote the two terms into brackets in this expression by 



d-{y) = S-(y)D 



(-y) 



S-(y)D-(y) 



(69) 



and 



d + (y) = S + (y)D 



(y)- 



(70) 



Taking the Fourier transform of d_ and using the explicit form of Pauli- Jordan 
distributions we obtain 



By virtue of the deltas and thetas, we have that q 2 = m 2 , q° < 0, (p — q) 2 = 
and q° —p° > 0. So, p = (p — q) + q is the sum of two 4- vectors, one timelike and 
the other lightlike, both in the backward light cone. So, p is timelike. Then we 
can choose a Lorentz frame such that p = (p°, ). We first consider the integral 
proportional to m in the integrand of this expression: 



For / this integral vanishes by symmetry considerations. For fi = the 
integral can be solved in a straightforward manner and results 



d_(p) 



j dqi(i + m)6 (-q°) S (q 2 - m 2 ) 9 (q° - p°) 6[{p - q) 2 ] . (71) 





The corresponding term in (71), in any reference frame, is given by 




(72) 



Now we calculate the term proportional to <j[ 2 in Q7ID: 
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Also this integral vanishes for \x 7^ v by symmetry reasons. By Lorentz invariance 
it must be proportional to g^ v . The proportionality coefficient is determined by 
saturation with g^ v . Thus 

V = \9r J dq q 2 e (-q°) 6 (q° - p°) 6 (g 2 - m 2 ) 5[(p - qf] . 

This integral can now be easily solved and gets the following contribution to 
integral flTT|): 

z 2 ' 
,2" 



rf-(p): 



'"" tfrPP ( 1 - ^ ) % 2 - m 2 )9 (-p°) 



(73) 



2(4vr) 3 - 

Combining (|72|) and (^) into (|71~D we get the final result for cL(p) = d-{p)\ + 
d_(p) 2 : 



d_(p) 



1 



:% 2 - m 2 )# ( 1 



(4tt) 3 ' 

Turning back to (|68D we have 

„2 



m 
p 2 



• (74) 



(4tt) 3 
fj(p), 



% 2 - m y (-/) ( 1 - ^ 



2m + p'[l + — 



(75) 



where we used the following algebraic relations for (3 matrices 



anc 



T We can show that this identity is valid using the usual scalar representation of section 2. As 
the r.h.s. is invariant under changes of representation, this result is general, i. e., representation 
independent. 
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The distribution d + (p) is calculated in the same way and results 



-e 2 (3»d + {p){3, 



6{p 2 -m 2 )6(p°) (l 



M (4tt) 3 
Now, (ff5|) and (|76|) together give 



m 
P 2 



2m + p* 1 



m 



■ (76) 



m 
P 2 



2m + /M 1 + 



m 



sgn(p°) . (77) 



We verify that the singular order of this distribution is u = 1. We then split 
it by using the central splitting formula (pIT| ) with uo = 1. For p timelike this yelds 



ri(p u ) 



/.+00 

2vr ^ ^ (fc°-20) 2 (p°- A;° + z0) 



(7* 



We turn (ffSD into this integral and solve it by standard methods 0]. For an 
arbitrary Lorentz frame we arrive at 



hip) 



le 



4(2tt) 4 

x < [log 1 1 - b 2 \ - iirsga (p°) 9 (p 2 - m 2 )] 



m I 1 



b 2 J 2 1 & 4 



± 

2b 2 



m — 



(79) 



with b 2 = ^2 • We can analitically extend this result for arbitrary complex p: 



le 



4(2tt) 4 



log(l-6 2 ) 



m 

2b 2 4 



Combining ( |79|) with (|75|), we obtain Tf eli (xi, x 2 ): 

Tj self (xi, x 2 ) = «e 2 : -?/>(xi)£ (x x - x 2 ) V> (^2) : , 

where 

g(p) = -z [ri(p) -rj(p)] 



JO) 



51) 



4(2tt) 



[log 1 1 — 6 2 1 -iitd (p 2 -m 2 )] 



m I 1 



b 2 



+ 



td- 



b 4 
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This result is free of infrared divergences, because we have used the central split- 
ting formula, which normalizes the solution at the point p = 0. Also, it is well 
defined on mass shell. Because u = 1 we must add to this result a first order 
polynimial in j> in order to get the general solution to the splitting problem. Thus, 

£(p) = ^{N|i-^ 2 |-^(p 2 -- 2 )][-(i-^)+^(i-^) 

-iL + Co + CiA- (83) 



2b 2 

We now consider radiative corrections due to the self energy insertions into 
the DKP propagator. The complete propagator is 

?c C o mp i(p) = T\p) + f %p)E (p) f%p) + f »£ (p) f »£ (p) T» + ... , 

where T c (p) is given by (|55|) . Formally summing this series we obtain 

compllPJ 1-E(p)f«(p) 

1 + m) — p 2 + m 2 



(27r) 2 m p 2 _ m 2 + - \j){J) + m ) - p 2 + m 2j s ( p ) 
We require that m be the physical mass of the DKP particle. This will be satisfied 
only if 

\J)(J> + m) — p 2 + m 2 ] E (p) = 0. 

p 2 =m 2 

Substituting ( ^3l) into this condition we obtain the following restriction on the 
arbitrary constants Cq and C\. 

C = m(l-Ci 



x 2 

This condition eliminates one of the arbitrary constants, say Co- The remaining 
one will be related by an Ward identity to another arbitrary constant that will 
appear in the vertex distribution. In the Appendix we derive this Ward identity 
by requiring gauge invariance of the theory. 
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Now we calculate the vertex distribution in the limit of zero tansferred mo- 
mentum by making use of the Ward identity d95|) deduced in the Appendix: 

1 d 



A"(p,p) 



(2tt) 2 dp, 



-E(p) . 



Substituting the explicit form of the self-energy ( [33l) into this identity we obtain 



ier 



4(2tt) 6 
2p" 



mb 4 
2pt* 



[log 1 1 - b' 



inO (p 2 — m 2 )] 



P.V. 



mb 2 
1 



6 2 

/3" 



1 



1 

2 V fr 4 



4) 



m 2 6 4 26 2 

This result is well defined at p — 0, but it is singular on mass shell p 2 = m 2 due 
to the logarithmic term. 

The above form of the vertex function suffices to study the physical meaning 
of the constant Ci, which will be done in conection to charge normalization. 
The physical charge is defined in the scattering of a scalar particle by an external 
electromagnetic field at low energies. Thus we must consider the contributions (in 
the limit of zero transferred momentum) to S matrix from the terms containing 
C\ in both self-energy and vertex distributions. Because of the above mentioned 
mass shell singularity, we must be care in taking the adiabatic limit. Making so, 
we can prove that all these contributions cancel themselves and conclude that this 
constant has no physical meaning. Nevertheless, it can be specified by requiring 
the vertex function to satisfy the central splitting condition, i. e. 



A> (0,0) = . 
Using this condition into ( |34D we obtain G\ = \. 



7 Concluding remarks 
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In this paper we have considered scalar QED based on Duffin-Kemmer-Petiau 
equation in the framework of Epstein-Glaser causal method. We have given 
the basis to construct the second order S matrix and calculate the lowest order 
distributions for Compton scattering, vacuum polarization and self-energy. By 
using the gauge invariance requirement we determined the vertex correction from 
the self-energy one. 

The starting point of the causal approach was the identification of the one- 
point distribution Ti(x), which was given by the interaction term in the La- 
grangian ([[3D of the theory, where all the fields entering that expression were free 
fields. Thus, the causal method specified completely the form of the interaction, 
giving us a non effective theory. After determining, by using gauge invariance, the 
finite normalization constant appearing in the two-point distribution for Comp- 
ton scattering, this distribution was identified with the propagator of the DKP 
scalar field. Then, in the causal approach we have recovered, in a natural and 
even simpler way, the basic quantities from which it is constructed the usual ef- 
fective theory. Namely, we identified iT c (x) as the effective propagator and — e : 
ifj(x)f3 IJ "ip(x) : (if) and if) being free fields) as the "vertex" (interaction). 

At one loop level we have calculated the scalar vacuum polarization tensor, the 
self-energy and the vertex correction. We have determined all physically mean- 
ingful finite normalization constants from physical requirements as symmetries 
and mass and charge normalization. Our results agree with that obtained in the 
context of the effective theory. To do a complete analysis at one loop level, we 
would have to calculate as well the singular order of the four point distribution 
for scattering of two scalar particles. This problem is fundamentally important to 
study the renormalizability of the theory and is under our investigation presently. 
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As future perspectives we can quote the use of the causal approach to study DKP 
field interacting with external gravitational fields. 
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Appendix: gauge invariance 

The gauge invariance requirement allows us to determine the remaining arbitrary 
constant C, in the two-point distribution for Compton scattering, and to find a 
relation (Ward identity) between the self-energy and vertex distributions, thus 
relating the corresponding arbitrary constants. 

As we saw, in the causal approach all the fields are free fields. Then, under a 
gauge transformation the electromagnetic field transforms as A^x) — > A^x) + 
d /Ji A(x), whereas the matter fields ip(x) and ip(x) remains unaffected. Here A(x) 
is a c-number scalar field that satisfies □A(rr) = (Lorentz gauge) and vanishes 
at infinity. 

The gauge invariance requirement amounts that in the adiabatic limit, g — ► 1, 
the S matrix must be invariant under such transformations, which implies the 
invariance of all the terms in the S matrix expansion (0). We write a generic 
n-point distribution normally ordered with respect to photon operators in the 
form 

n 

T n (xi, x n ) = y ] ^] tki...ki (^i' Xn ) • A«i ( x ki)---^m ( x ki) '■ j 

1=0 l<ki<...<ki<n 

where x n ) contains only scalar operators and it is the sum of all 
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graphs of order n with only I external photon lines, attached at the vertices 
, ■ ■ ■ , Xk r The external scalar operators are arbitrary. Applying the gauge 
invariance requirement we arrive at the following condition | |2T|j : 

^C;.X(x 1 ,...,x„) = 0, (85) 

kj 

for all 1 < I < n, all 1 < j < I, all 1 < k\ < ... < k t < n and all (x 1 , x n ) G R 4n . 

When we apply this condition for the two-point Compton distribution, we get 
at the condition 

dlQ^( Xl ,x 2 )=0, (86) 

where 

Q fM/ (x 1 ,x 2 ) = ■.^(x 1 )/3H l (x 1 ,x 2 )/3^(x 2 ): 

+ :^(x 2 )p u t n (x 1} x 2 )P^(x 1 ) : 
= Q Ufl (x 2 ,xi) . (87) 

and means derivative with respect to x±. 

Turning (|8~7D , ([49]) and (|5l| ) into this last equation, the l.h.s. gives 



dlQ^{ Xl , x 2 ) = - : d^(x 1 )^[S F {x 1 - x 2 ) - CS{ Xl - x 2 )](3 v il>{x 2 ) : 

- : ^(x 1 )^[d,S F (x 1 - x 2 ) - Cd,5( Xl - x 2 )}(3^(x 2 ) : 

- : ^(x 2 )(3 u [-d f ,S F (x 2 - xi) + Cd fl 5(x 2 - x 1 )\^{x 1 ) : 
-■.^{x 2 )(3 v [S F {x 2 -x x )-C5{x 2 -x 1 )\(3»d^{x x ): . (88) 

From the equation of scalar Feynman propagator 

(□ + m 2 )A F (:r) = 5(x) , 

the relation (f|) and definition (|C]), we have 

Fi x) = mS F (x) + -38(x) . (89) 
m 
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Using this into (88), condition (36) yields 



C 



I 

m 



(90) 



where / is the 5x5 identity matrix. 

Now we consider the three-point distributions corresponding to the vertex 
correction. We must consider four graphs, that are analogous to that we find in 
the vertex correction of spinor QED. Considering the photon line attached at the 



vertex x 3 , the distribution t 3 (xx,x 2 ,x 3 ) we must enter into condition fl85|) is 



t 3 l (x 1 ,x 2 ,x 3 ) 



- : ip(x 1 )A (i (x 1 - x 3 , x 2 - x 3 )ip(x 2 ) : 

+ : #£!)£(£! - x 2 )T c (x 2 - x 3 )^{x 3 ) : 

+ : ^{x 3 )^T c {x 3 - Xi)E(a; 1 - x 2 )^{x 2 ) : 

+ : : g pu D^ Xl - x 2 )\T»(x 2 - x 3 ) 

+Xi < — > x 2 , 



where T c is the "effective DKP propagator" given in (p3|) . Thus we have 

d 3 ^(x 1 , x 2 , x 3 ) = : ij}{xx) [-9jA M (xi - x 3 , x 2 - x 3 )] ip(x 2 ) : 

+ : ^(xOE(x! - x 2 )dl [T c (x 2 - x 3 )f3^{x 3 )} : 
+ : &l $(x 3 )^T c (x 3 - x x )] - x 2 )ij(x 2 ) : 
+ : ^(xO : g pv Dl (x x - x 2 ) [dfil^(x 2 - x 3 

+Xi < — ► x 2 . 



For the brackets in the second and third lines we have 



(91) 



dl [T c {x 2 - x 3 ) n(x 3 )] 



x 3 )^T\x 3 - xx 



i5(x 2 - x 3 )i[)(x 3 ) ; 
-itp(x 3 )5(x 3 - xx) , 



where we have used the DKP equation and the fact that T c is the corresponding 
DKP Green function (see eq. fl5~3])). Turning these results into fl9~ll) and taking 



32 

into account the delta distributions, we obtain 

d$t%(xi,x 2 ,x 3 ) = : 4>(xi) {-^A M (xi - x 3 ,x 2 - x 3 ) 

+i5(x 2 - x 3 )S(xi - x 2 ) - i5(xi - x 3 )E(xi - x 2 )} ip(x 2 ) : 
= . (92) 

This condition will be satisfied only if the term into curl brackets vanish. Thus 
we arrive to the desired Ward identity: 

-df l A fl (x l -x 3 ,x 2 -x 3 ) + i5(x 2 -x 3 )Yj(x l -x 2 )-i5(xi-x 3 )Yj(x 1 -x 2 ) = . (93) 

Denoting 



we can write 



yi= xi- x 3 and y 2 = x 2 - x 3 , 



^A" (y 1 ,y 2 ) = -(d 1 + d 2 )^(y 1 ,y 2 ) 



dx 3 



and equation ( |9"3] ) becomes 

(d 1 + d 2 )^ (y h y 2 ) + i5(y 2 )E( yi - y 2 ) - i5( yi )E( yi - y 2 ) = . 
Taking the Fourier transform^ of this equation we arrive at 

(pi + p 2 ) m a> - {s(pO - s( P2 )} = o . 

Introducing the notation q = —p 2 , p — Pi and redefining as 
A^(p, -g) - A"(p, g), 

■^Our convention to the Fourier transform of a distribution F(yi,y2) is 
%i,p 2 ) = ^ 1 <%*/2 e^+P^F( yi ,y 2 ) 
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we arrive at 

(p-q)^(p,q) = -^ {%)-%)} . (94) 
Taking the limit p — > g we have 

x " (p - pH (^l: £(p) - (95) 

which is the form of the Ward identity we use to determine the vertex correction 
in Section 7. 
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